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Abstract 

A linear bubble model of grain growth is introduced to study the condi- 
tions under which an isolated grain can grow to a size much larger than the 
surrounding matrix average (abnormal growth). We first consider the case of 
bubbles of two different types such that the permeability of links joining unlike 
bubbles is larger than that of like bubbles (a simple model of grain boundary 
anisotropy) . Stable abnormal growth is found both by mean field analysis and 
direct numerical solution. We next study the role of grain boundary pinning 
(e.g., due to impurities or precipitate phases) by introducing a linear bubble 
model that includes lower and upper thresholds in the driving force for bubble 
growth. The link permeability is assumed finite for driving forces above the 
upper threshold, zero below the lower threshold, and hysteretic in between. 
Abnormal growth is also observed in this case. 
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I. INTRODUCTION 



We use linear bubble models of grain growth, originally developed to study self-similar 
particle coarsening or the development of texture, to investigate possible causes of abnormal 
grain growth. While the bubble model is a simple idealization of grain growth, it has the 
advantage that the results obtained are not limited by the mean field approximation inherent 
to other existing treatments of abnormal growth. 

In normal grain growth thermal annealing of a polycrystalline material results in self- 
similar coarsening driven by excess free energy reduction. An invariant distribution of scaled 
grain sizes develops, with an average grain size that grows as a power law of time with a 
characteristic exponent of 1/2 In abnormal grain growth, on the other hand, a few 

grains grow to a very large size relative to the average matrix. In some cases a stable 
operating state is achieved characterized by a constant ratio of abnormal grain sizes to 
matrix average. 

Different mechanisms have been proposed for abnormal grain growth, and some tested 
against Monte Carlo simulations of discrete, lattice models of a polycrystalline material. 
Defect induced strains can induce isolated grain growth p, as well as the same capillary 
forces responsible for coarsening when anisotropy of grain boundary energies or mobilities 
exist . The conditions for abnormal growth due to variable surface energies or mobilities 
were recently examined in ref. within a mean field treatment of the matrix grains. For 
the case of a single grain with boundary properties that differ from those of the surrounding 
matrix, it was found that a higher boundary mobility generally promotes abnormal growth 
whereas a higher boundary energy constrains it. The detailed behavior can be quite complex 
depending on the ranges of the model parameters chosen. It includes abnormal growth only 
up to a limiting grain size, or lower bounds in the initial size of the grain for abnormal 
growth to occur. 

Abnormal grain growth has also been shown to occur when grain boundaries pin due 
to, for example, existing precipitate phases or other defects. Simplified models have been 
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proposed that introduce grain boundary drag forces that lead to ultimate pinning (Zener 
pinning) |]rU|,|TT|, while the role of thermal fluctuations to overcome pinning has been analyzed 



by Monte Carlo simulation 



Linear bubble models of grain growth were originally introduced by Hunderi et al. [|13 
Their results showed self-similar coarsening with parabolic growth kinetics for the average 
bubble size. The effects of grain orientation and anisotropic boundary properties have also 



been introduced into these models by Novikov , and later by Abbruzzese and Liicke |T5 
to study the development of texture. Both in Novikov's work and in later work [p!^ -[T8| two 
types of bubbles A and B were considered to represent the idealized situation of only two 
different grain orientations. The mobility of unequal A-B boundaries was assumed to be 
larger than the mobility of either A-A or B-B boundaries. The conditions under which a 



steady state distribution is reached in this binary case was studied in ref. . 

We consider in this paper a linear bubble model of grain growth to investigate the 
conditions that could promote abnormal growth in the model. Section |l| presents a mean 
field analysis along the lines of the mean field treatment of grain growth by Rollett and 
Mullins 0. We consider a linear chain of bubbles of two types to model grain boundary 
anisotropy and, consistent with their analysis, show that abnormal bubble growth is possible 
when the permeability of links separating unlike bubbles is larger than that of like bubbles. 
This mean field analysis is complemented by a direct numerical solution of the model that 
confirms the mean field predictions regarding abnormal bubble growth: parabolic kinetics 
for both the abnormal bubbles and the matrix, and a constant value of the ultimate ratio 
of abnormal to matrix bubble radii. 

We then explore in Section |T| a different mechanism that can lead to abnormal growth 
even for the case of a single bubble type (i.e., in the absence of mobility anisotropy). We 
model grain boundary pinning by introducing a finite threshold in the driving force for 
bubble growth. Links between adjacent bubbles are open if the driving force exceeds an 
upper threshold, and closed if it falls below a lower threshold. In between the two thresholds 
we assume hysteretic behavior of the link permeability. A numerical solution of the model 



shows that abnormal growth follows from an initial steady state distribution of bubble sizes 
depending on the values of the upper and lower thresholds. While the upper threshold largely 
determines the subset of bubbles that can grow against the matrix, we show that there is 
a sharp transition in behavior depending on the value of the lower threshold. Abnormal 
growth is observed below a critical value, with the average radius of the bubbles growing 
linearly with time. Above this value, the bubble distribution freezes after an initial transient 
(growth stops). 

II. ASYMMETRIC LINEAR BUBBLE MODEL 

A linear bubble model of abnormal grain growth is introduced to address the relationship 
between stable abnormal growth and anisotropic grain boundary mobility. The analysis is 
motivated by recent research that involved the idealized situation in which a single isolated 
grain A grows in a matrix of B grains . Under the assumptions that the A — B boundary 
has a different energy and mobility than B — B boundaries, that the boundary vertices are 
in equilibrium, and a mean field treatment of the B grain matrix, it was concluded that 
abnormal grain growth is to be observed when unequal boundaries have higher mobility 
than equal boundaries, whereas higher surface energy of unequal boundaries relative to equal 
boundaries constrains it. The linear bubble model described here allows us to extend these 
results beyond the mean field approximation for the matrix. We find that the conclusion that 
abnormal growth occurs when unequal boundaries have higher mobility also holds in this 
case, and that a mean field prediction of the ultimate size ratio is in reasonable agreement 
with the results of the numerical calculations. 

We consider a set of spherical bubbles of radii Ri, i = 1, . . . , N, forming a linear chain 
with periodic boundary conditions. The temporal evolution of the linear bubble model is 
defined by the following set of equations. 




(1) 
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where Mij is a permeability coefficient between bubbles i and j (the analog of the mobility 
in the grain boundary case). We ffist consider in this section the case of two types of grains, 
A and B, such that Maa = Mbb = 1? and define /i = Mab/Mbb- A general property of 
Eq. dlD is the existence of a conserved quantity, J2iLi Ri, which is independent of time. 

Numerical results for the symmetric case fi = 1 were given in ref. For a random 

initial distribution of bubble radii, the ensemble coarsens through growth of bubbles larger 
than a time dependent critical radius, and shrinkage and disappearance otherwise. Following 
an initial transient, the configuration reaches a stationary self-similar state. In it, consecutive 
configurations of the coarsening structure are geometrically similar in a statistical sense. As 
a consequence, any linear scale of the structure (i.e., the average bubble radius) grows as a 
power law of time 



where ( ) denotes the configuration average, and to is some time in the self-similar regime. 
Before presenting the results of our numerical calculations for /i > 1, we discuss a mean 



show that similar conclusions follow. We then obtain a numerical solution of Eq. ([^), and 
demonstrate that, in agreement with the mean field results, the linear bubble model does 
lead to abnormal grain growth when /i > 1. We also show that the ultimate size ratio 
between the abnormal grains and the matrix is quite close to that predicted by the mean 
field analysis. 

Consider a bubble of type A in a long chain of B bubbles, and that the AB links have 
a mobility Mab 7^ Mbb, with the mobility ratio /i = Mab/Mbb- We calculate the time 
dependence of = Ra/ (Rb) by using a mean field approximation to the evolution of the B 
bubbles. We start from. 



{Rit))'-{R{to))' = Cit-to), 



field treatment of the linear bubble model (Eq. (|1])) along the same lines of ref. 0, and 





Since the A bubble has two B bubbles as neighbors, one has. 
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{Ra\Ra) = 2AitB 



Rb Ra 



2M 



AB 



a 



(3) 



_{Rb) Ra 

where the second equahty follows from assuming self-similarity of the matrix bubbles and 
a = {Rb){1/Rb) = {Rb)/Rc ^ 1.1927 (cf. Appendix 0). The mean field treatment of the 
B bubbles outlined in the Appendix gives for the critical radius Rc = Mbb/'^Ro and hence 
from the definition of a we find, 

d{RB) MBBd^ 



dt 2{Rb) ■ 
Substituting Eqs. (||) and (§) into Eq. (|^) gives, 



(4) 



where, 



i^l^) = T^\^G{fi,uj), (5) 
{Rb}^ 



G(/.,.;) = 2/.(«-i)-^. (6) 



This latter function determines the sign of {uj\uj), and therefore whether the A bubble grows 
or shrinks relative to the coarsening B matrix. 

For fi = 1 the function G is everywhere non-positive. For fi > 1 there is a range of values 
of u for which G is positive, and in particular a stable fixed point at some u = that 
corresponds to steady abnormal growth. Figure |l] shows the phase space plot of u for a = 1 
and fi = 1.5. For values of uj from roughly 1 to 4, cj > so that a bubble of type A in this 
range would grow relative to the matrix of B bubbles. However, if the ratio u) exceeds 4, the 
larger bubble would shrink back to the fixed point. This root of G is a stable fixed point. 
The other root ~ 1 is not stable. The upper root of uj) is given by, 
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w+ = - -/ij . (7) 

The range of relative growth is given by the difference between the upper and lower roots 
of Eq. i), 

= -^fi'-fi. (8) 
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We next compare the results of the mean field calculation to a direct numerical solution 
of the set of equations (|I]). We only describe the algorithm briefiy, further details can be 
found in ref. |T^. We consider a large number of bubbles N = 2 x 10^, and impose periodic 



boundary conditions such that RN+i{t) = Ri(t). We initially place 20 equally spaced bubbles 
of type A in a matrix of — 20 bubbles of type B. The initial sizes of B bubbles are 
distributed according to the mean field distribution (Eq. (|A9|) with {RB)(t = 0) = 5). The 
initial radius of the A bubbles is fixed at RA(t = 0) = 3{RB){t = 0) for the results presented 
in this section. A wide range of initial ratios has been investigated with identical results. 
We also set Mbb = 1. A lower size cut-off Rmin is introduced for numerical reasons so that 
any bubble for which Ri{t) < Rmin during the course of the calculation is removed, and the 
two adjacent bubbles redefined as neighbors. The value of Rmin = 0.28 is chosen so that no 
bubble can shrink to zero in At = 0.02, the time discretization used to integrate the system 
of equations (|T|). The averages shown refer only to averages over the configuration. We have 
not performed additional averages over independent initial conditions as the large number 
of bubbles considered appears to be sufficient for the required statistical accuracy. 

Figure || shows our results for (/^^(t)) as well as {Ra), where the latter is an average 
over the 20 bubbles of type A. The average radius of the matrix bubbles (/^^(t)) exhibits 
power law growth with an exponent of 1/2, in agreement with the mean field prediction 
of Appendix |^. The figure also shows two least square fits to obtain the corresponding 
amplitudes of the power laws which are used to calculate the ratio {Ra) / (Rb) at long times. 
Figure ^ shows our numerical results for the ratio of amplitudes for a range of mobility ratios 
/i, and compares them with the mean field prediction given by Eq. (|^. 

In summary, a mean field treatment of the linear bubble model with unequal boundary 
mobilities predicts that abnormal bubble growth will occur for /i > 1 with an ultimate size 
ratio of ujj^. The numerical results confirm power law growth in time of both A and B average 
radii, with an exponent of 1/2. The numerical results for the ultimate size ratio u;+ are also 
in excellent agreement with the mean field prediction. Clearly, bubble size correlations that 
are not taken into account in the mean field treatment must only introduce small corrections. 



III. SYMMETRIC CASE WITH A MOBILITY THRESHOLD 



We investigate in this section a different mechanism leading to abnormal grain growth 
even in the absence of any mobility anisotropy. We hypothesize that if a finite threshold to 
grain boundary motion exists, then it is possible that a large fraction of the matrix grains 
would remain immobile, except for those that were sufficiently larger than their neighbors 
so that the local driving force for growth exceeds the given threshold. The excess energy 
that is contained in the initial particle distribution would then be relieved mostly through 
size increases of the larger grains at the expense of a largely immobile, high energy, matrix 
distribution. 

In order to investigate this possibility within the linear bubble model introduced in 
Section we consider an ensemble of like bubbles and introduce two threshold values for 
the mobility M in Eq. (|lD. Let Ap = 1/Ri — be the local driving force associated 

with the i — th link, and Api < Ap„ the low and high driving force thresholds respectively. 
We define M = 1 if |Ap| > Apu, and M = if |Ap| < Api. We also assume a hysteresis 
loop in Api < \Ap\ < Apu with M = in the lower branch and M = 1 in the upper branch. 
Therefore a link remains closed (M = 0) until \Ap\ across the link exceeds Ap^. Once the 
link is open (M = 1) it remains open until |Ap| falls below Api. Finally, when a bubble 
radius falls below Rmin, so that the bubble is removed from the distribution, a new link 
between the new neighboring bubbles is made, and its mobility is assigned to be 1 unless 
\Ap\ < Api. 

We have used the same numerical algorithm described in Section H to integrate the 
system of equations (|l|) with the mobility thresholds just introduced. In this Section we 
consider an ensemble of = 10® identical bubbles, initially distributed according to the 
mean field result, Eq. ([X9|), with {R{t = 0)) = 5. Although all bubbles are identical in 
the present case, it is convenient for the sake of the discussion to refer to those that grow 
relative to the average as A bubbles, and as B or matrix bubbles to the rest. 

Figure § shows our results for the bubble radius distribution function for a representative 
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set of parameters Api = 0.05 and Apu = 3.0. Even though all bubbles are identical and 
follow the mean field distribution at t = 0, the largest bubbles in the initial ensemble grow 
while most of the rest remain stagnant. This figure shows the radius distribution p{R) (with 
the main peak near R = {R) suppressed for clarity) starting at t = 1000 all the way up to 
t = 25000 in increments of 2000 time units. It is clear from the figure that a small subset of 
the initial distribution grows as indicated by the successive peaks of p{R) at large R. 

Not all possible combinations of Api and Ap^ result in abnormal growth however. First, 
there is an obvious upper bound for Ap^ given the initial radius distribution, and it corre- 
sponds to the driving force Ap between the largest possible bubble Rmax — 8.4 for our initial 
distribution, and Rmin - 0.28. We find Ap^ax = ^/Rmin - ^/Rmax - 3.42. If Ap„ > Ap^ax 
no bubble will grow. 

For fixed Apu < Apmax, a bubble will grow (call it A) at the expense of a B neighbor 
when Ra > {1/Rb — ApuY^ ■ Therefore for a given initial distribution the value of Ap^ 
determines the range of radii of bubbles expected to grow. Once a given A bubble starts 
growing, it will only stop if it encounters a bubble B such that 1/Rb — ^/Ra < ^Pi- 
Equivalently, whenever an A bubble encounters a bubble of radius Rb = {\/ Ra + Api)^^ or 
larger, growth will stop. If one further assumes that growth of A has already occurred for 
some time so that Ra is sufficiently large, then this condition is approximately Rb ~ 1/ Api 
independent of Ra-, relation that can be used to defined a critical value for growth Ap^. If 
Api > Api there is a nonzero probability that a large and growing A bubble will become the 
neighbor of a B bubble that is sufficiently large to stop growth of the A bubble. If the matrix 
has remained approximately stagnant, this critical value can be obtained from Rmax — 8.4 
as Api = 1/Rmax — 0.12 (in practice, the numerically sampled initial condition typically 
has Rmax — 7.9 or Apf ~ 0.127). This behavior is observed numerically and is illustrated in 
Fig. 1^. The figure shows the radius of the largest bubble in the ensemble as a function of 
time. The upper cut-off in all the cases shown is fixed Apu = 2.0, and the figure shows the 
results for a range of values of Api. The value Api ~ 0.130 marks the transition between 
abnormal growth and an ultimately frozen configuration. Note that the transition is quite 
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sharp as a function of Api. Identical numerical results concerning this transition as well as 
the same critical value Ap^ have been obtained for Ap^ = 2.5 and Ap^ = 3.0. 

In summary, the upper cut-off Ap^ determines the fraction of the ensemble that can 
grow, and therefore the degree of stagnation of the matrix. Once abnormal grain growth 
has started, the value of the lower cut-off Api (and the amount of growth in the matrix, if 
any) determines whether abnormal growth continues or rather the system reaches a frozen 
configuration. 

We finally mention that while abnormal growth occurs, the typical radius Ra of the large 
particles is expected to grow linearly with time. In mean field, a given A bubble will have 
two B bubbles as nearest neighbors and therefore 



While bubble A grows successive B neighbors will shrink to zero and be eliminated from 
the ensemble. Therefore the growth of A can be estimated by averaging Eq. (|^) over the 
distribution of B, and when the matrix is almost stagnant, over the initial distribution of 
bubble radii. In either case, {1/Rb), where ( ) denotes average over the configuration, will 
be constant (or changing very slowly compared with the rate of growth of the A bubble), 
so that for sufficiently long times {1/Ra ^ {^/Rb)) dRA{t)/dt is approximately constant. 
This is the dependence shown in Fig. (|^) in the cases in which abnormal grain growth is 
observed. The value of the slope of the line RA{t) versus t can be straightforwardly obtained 
by computing the average {Mab/ Rb) over the initial distribution of matrix grains. 
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APPENDIX A: MEAN FIELD CALCULATION FOR A LINEAR BUBBLE 

MODEL 



Consider a set of spherical bubbles arranged along a line that model a set of neighboring 
(identical) grain boundaries. The equation of motion for their radii is, 

^ - M + — —] (Al) 

dt \Ri+i Ri~i Ri ) 

The mean field approximation for the average (ignoring correlations) is, 

(i?,|i?,) = 2M(i---i), (A2) 

where one defines a critical radius through 1/Rc = (l/R). Define now a reduced radius 
r = R/Rc. Equation (^) can be written as, 

Rc{Ri\Ri) =2M(^l-^y (A3) 



As is standard in the analysis of steady state solutions for the averages |T9|, one first defines 
the quantity, 

V — RnRr — ~ — ; — , 

^ 2 dt 

so that Eq. (^3) can be written as. 



Rl{r\r) = f{r)-ry, (A4) 

where we have defined /(r) = 2M(1 — 1/r). This equation is a particular case of Eq. (7) in 
ref. [^. The nodal curve defined by (r|r) = is thus given in our case by. 



y = 2M(--l). (A5) 



According to the classical mean field treatment of Lifshitz and Slyozov [^], there exists a 



stable operating point of the reduced particle size distribution determined by Eq. (|A1|) that 
corresponds to the maximum of the nodal curve y = Umi so that the distribution of reduced 
radii r extends from r = to a sharp cutoff r = r^- For our particular form of the nodal 
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curve, Eq. (|X5|), we have = 2, and ym = M/2. A statistical self-similar distribution is 
reached with this value of y, and from its definition, we have, 

dRc _ M 

that after integration leads to the asymptotic parabolic growth law, 

Rl{t)-Rl{t,) = ^{t-t,), 

where to is some time within the self-similar regime. 

The distribution of reduced particle sizes can also be computed by using our result for 



/(r), and Eq. (14) in ||T9[. Define the function F{r) = t(r|r), which satisfies in the steady 
state, 

F{r) = = _(^. (A6) 

The general solution of the continuity equation for n{r,t), the number particle density, is 
given by |T9[, 



nir,t) = ^^ [t-Qir) 



where is an arbitrary function and where. 



/•r fjrf' 

JQ F[r) 

Substitution of Eq. (^) into Eq. (^) yields, 

e(r)/2 = -ln(2 - r) - + ln2 + 1. (A8) 

With this result, the normalized probability distribution function P(r) = n{r,t)/ J n{r,t)dr 
is time independent and given by (Eq. (20) in ref. [0), 

P(r) = ^^^^e-2/(2-'-), 0<r<r^. (A9) 

We note that the upper cut-off is = 2, that the maximum of P(r) occurs at r = ^/2, and 
that the average reduced radius is given by (r) = {R) /R^ ~ 1.1927. 
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FIG. 1. Phase plot uj versus uj where to = Ra/{Rb) is the ratio between the radius of the A 
and matrix bubbles. G > corresponds to ratio growth, and G < otherwise. The plot shows 
two fixed points {G = 0) at two different values of w. The smallest of the two is unstable, and the 
largest, denoted by u;+ is stable. This is the expected operating point of the model and corresponds 
to a fixed size ratio between Ra and {Rb)- 
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FIG. 2. Average radius of large bubbles {Ra), and average radius of the matrix {Rb) as a 
function of time for a mobility ratio /i = 1.5. Both A and B bubbles exhibit average parabolic 
growth to a very good approximation as shown by the fits (solid lines). The amplitudes of the 
term t^^"^ are used to estimate the quantity uj-^ shown in Fig. |3| for each value of ^. 
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FIG. 3. (o), values of uj-^- obtained by numerical integration. For each value of the mobility ratio 
fi the model equations are integrated in time, and the value of a;+ estimated by fitting parabolas 
to the average radii as shown in Fig. |2|. The solid line is the mean field prediction, Eq. (^. 
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FIG. 4. Probability distribution function of particle radius as a function of R for several times 
ranging from t = 1000 to t = 25000 in increments of 2000 time units. The main peak of the 
distribution has been removed for clarity. Each of the peaks shown at large R corresponds to a 
specific time, and they are ordered from left to right according to increasing times. The distributions 
show the existence of a small set of large bubbles that grow to a size much larger than the matrix 
average. 
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FIG. 5. Largest bubble radius in the distribution as a function of time. The upper threshold 
Apu = 2.0 in all the cases shown, and the values of the lower thresholds are indicated in the 
figure. The three largest values of Api = 0.131, 0.150 and 0.200 lead to an asymptotically stagnant 
configuration (left axis), whereas the values Api = 0.100,0.125 and 0.130 lead to abnormal growth 
(right axis). 
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